“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1969 


Optimal design of trussed structures. 


Crockett, William Arthur 


Naval Postgraduate School 
http://ndl.handle.net/10945/12774 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 
1969 
CROCKETT, W. 


OPTIMAL DESIGN OF TRUSSED STRUCTURES. 
by 


William Arthur Crockett 











DUDLEY KNOX LIBRARY 


NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CA 93943-5104 


United States 
Naval Postgraduate School 





THESIS 


OPTIMAL DESIGN OF TRUSSED STRUCTURES 


William Arthur Crockett 





October 1969 


This document has been approved for public re- 
Lease and sale; its distribution 4 unlimcted. 








Optimal Design of Trussed Structures 


by 


William Arthur_Crockett 
Major, United States Army 
B.S., University of Texas, 1960 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN OPERATIONS RESEARCH 


from the 


NAVAL POSTGRADUATE SCHOOL 
Octoscr F169 


ABSTRACT 
Structural design problems can be considered to be optimization 

problems because a design is sought which is optimal by some criterion 
subject to limitations on size, behavior, or other aspects of the struc- 
ture. Under certain conditions, such problems may be solved by con- 
ventional mathematical programming techniques. The minimum weight 
design of a statically indeterminate three-bar truss is used to illustrate 
the application of the "Sequential unconstrained minimization technique" 
of Fiacco and McCormick to the optimal design of trussed structures. A 


suggested computational procedure and computer program are included. 
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Lee INTRODUCTION 

The desired conclusion of the structural design process is a structur- 
al configuration which performs its intended function efficiently. The 
method by which the final design is determined usually includes the 
establishment of an initial configuration from which the final design is 
evolved through a process of analysis and redesign. 

The electronic computer has permitted an acceleration of the design 
process by enabling the rapid analysis of the mathematical models assoc- 
iated with the structures under consideration. In order to more efficiently 
use this capability L. Schmit [1] has suggested a rational approach to the 
design process called "systematic structural synthesis." This approach 
involves the systematic evaluation and modification of a large number of 
trial designs until the optimal design is obtained. 

The basic problem is to design a structure which is optimal by some 
criterion subject to a set of requirements which specify acceptable limits 
on the behavior, size, or other aspects of the structure. Such problems 
are known to the operations analyst as mathematical programming problems, 
the general form of which is as follows. Find a vector X which 


minimizes F(X), 
suject to (Ha(x 
G, (X 


) =)0p i= lo2 aes 
|e iO ah 2 ae x OYA 


ali a « e 
Here, X = (xX,,Xe,...,X_) is an n-dimensional column vector; 
n 


F(X) is the objective function; and the relations H,(X) = 0 and G,(X) 20 


represent the constraints on the problem. 


If F(X) is convex in X and the constraint region is convex, i.e. the 
set of alll solutions tojb edges. ..,m and G,(X)=-0, j =1,...7p 
forms a convex set, the problem is called a convex programming problem. 
If a problem can be classified as a convex programming problem then 
there are solution procedures which are guaranteed to find the global 
optimum. 

Several techniques which have been devised to solve this special 
case of the general programming problem are Rosen's "gradient projection 
method," Zoutendijk's "method of feasible directions," and Kelly's 
"cutting-plane" method. These procedures are summarized by Hadley in 
Ref. 7. Another convex programming technique called the "method of 
alternate steps" was suggested by Schmit [1] as a means for converging 
on the optimal structural design. 

A recently developed convex programming technique which appears 
to be especially promising is the "Sequential unconstrained minimization 
technique for convex programming with equality constraints" of Fiacco 
and McCormick [2]. By this method, the objective function and con- 
straints are dealt with simultaneously by combining them into a single 
function. This function is minimized for different values of an auxiliary 
parameter thereby generating a sequence of feasible solutions that con- 
verge to the optimal solution. 

If the convexity restrictions for the constraint region or the objective 
function are relaxed a convex programming solution technique may still be 


used but there is no guarantee that a solution to the problem is other than 
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a local minimum. For most practical problems, it is no less valuable to 
obtain information about local optima in the absence of a global solution. 
In general, structural design problems do not meet the convexity 
restrictions for convex programming problems. One such problem is the 
determination of the cross-sectional areas of the members of a statically 
indeterminant, coplanar, three-bar truss such that the total weight of the 

truss is minimized subject to constraints on stress and displacement. 
This structure is illustrated in Figure 1. Procedures which can be used 


to solve this problem can also be applied to other trussed structures. 


leh 


IT. OBJECTIVE AND SCOPE 


The purpose of this study was to illustrate the application of the 
Fiacco-McCormick technique to the optimal design of trussed structures. 
Since minimum weight has usually been the major goal in the development 
of aircraft structural design techniques, it was used as the criterion for 
optimality. 

One example of the determination of the minimum weight design of a 
three-bar truss is included with results and the computer program. Also 
included is a brief discussion and outline of the solution procedure used 
to solve the problem. This procedure can only be applied to problems for 


which constraints can be formulated in completely analytical form. 


RZ 


Ill, THE THREE-BAR TRUSS PROBLEM 
The three-bar truss problem of Schmit (l]Jis illustrated in the following 


sketch. y 


i 


4 


u 





xj 


Figure 1. The Three-Bar Truss 


The function of the truss is to transmit the load P, from the point 
of application at joint A to the support B. The deflections of the joint 
Ain the x and y directions due to the load P, are Le and ae respectively. 
For a coplanar system of forces, one of the conditions for equili- 


brium of a rigid body in one plane is that the algebraic sum of the 


ke 


components of all forces in each of two mutually perpendicular directions 
in the plane of the forces must separately vanish. Therefore, 
LF = p,, cosf, + pe, cosBz, + pg, cosB3s + P, cosa, = 0 
and 
xP PP) th Pay ciies 1 Day oll b siiice 80 
where p,, is the force component in the direction of the ith structural 
member due to the jth load E , 
The stress in the ith member due to the jth force is defined as 
O15 = Pyy SA 
where Ais tilegcrocs-scclonal ared of the gen member sa lm@ered Uili- 
brium equations can then be written as 
A,0,, cosB, + Ag,,; cosBz + Ag 03, cosh; + P, cosa, = 0 (1) 
and 
A, 0,, sinB, + Ag@g, sinBy + AgOg, sinB, - P, sina, = 0 (2) 
Assuming small displacements, the change in length of the ith 
structural member due to the jth load is 
6.,=-uU_., COSsR=nI eesimee. (3) 
1) xj B, vi { 
A relationship between stess, displacement, and temperature change can 


be established if is noted that re may also be expressed as 


OR = 
eee 
Oy = A, E, +a, 47T,1,. (4) 


where jis the length of the ith member, a, is the mean cociiicicnives 
thermal expansion of the ith member, E, is the modulus of elasticity of 
the ithememiber, andgAT, is the temperature change imposed externally on 


the ith member. 


If the temperature is assumed to be constant, equations (3) and (4) 
7, 


can be combined to yield 


Pil, 
ae spb, + Simp, = 70 

Aye, Des cosh, ty By 

and since 0,, = p,,/A, and 1, = N/sing&, , this equation becomes 


Orr 
<i n 7 = 2 
Esme) ae i cosh, + u sinf, = 0. (S) 


For purposes of this illustration, it is assumed that all geometric 
parameters of the truss are given except for the cross-sectional areas, 
It is further assumed that in addition to non-negativity restrictions for 
the A, values, the behavioral limits for the structure are specified by 
upper and lower bounds on the stresses and displacements. Then, for 
the jth load condition, 

liggg= 84, =Uyy . b= Ds ee 


a ie ee 
XJ XJ XJ 


and 


lees tl See ee 
y) y) y) 


If o, denotes the specified density of the material in the ith struc- 
tural member then the weight, W, of the truss is 
W =A, 1, 0, + Az le 02 + Ag ls ps 
Since 1, = N/sinf, , this equation can be rewritten as 


“= A N py 4 Az N Oo 4 As N fs 


sing, sins sink, 


which is linear in A and is therefore convex, 


Ike: 


The mathematical programming problem is then to fine A, , Ap, As, 


the cross-sectional areas of the structural members which minimize 








N* py N* be N*03 
= oe 58 
iad Ay sinf, 2 sings As sinfs 


subject to the conStraints 
Ayo, cosB, + Az, cosh, + A303, cosh, + P, cosa, = 0 
A, O,, sinB, + AzO2, sinbz + A303, sinBs - P, sina, = 0 


N 


——e a S65, + S] = 
On, ee Hee By bo inf, 0 
N 
Gg Se 7 Ue eCos +u .sin = 0 
eee leyeeuijoy ey xj Pe yj Bz 
N 
i. COS. + UU. Sti = 0 
03; ESinBs xj bs yj Bs 
and 
Le = On s U;, Peas E c = u 4 < U 4 A < u ‘ < U A J A, 20 
x) x) x) Yj y) yj 


fori=1,2,3 andj =1,2,...,k where k denotes the number of distinct 
loading conditions. It is assumed that the loading conditions are not 
applied simultaneously. Note that inequalities such as L,, < 0,, < Uj, 
can be rewritten as 

Oe ee = a0, 

poem wera = AUP 

Now, in this problem statement, the variables to be determined are 

actually A, ,Az 1Ag.Oq5, i-l 2,70, j=l, 2 kh one Cee he , J =e, .. ole 
It is therefore apparent that the first two equality constraints are non- 
linear because of the cross products 0,, A, and hence the constraint 


region iS noneconvex., 
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IV, SOLUTION PROCEDURE 
Recall that the general form of the mathematical programming problem 
is tO determine the vector <~ wien 
minimizes F(X), 


1,2,...,M; 


I 


subject to H, (X) =0, i 


aia) =e 
where a =0@  Xetjiuum ee is an n-dimensional column vector. 


The basis of the Fiacco and McCormick technique is the P-function 


which is defined as follows for this general form. 


p 
x OS 


M 3 
4+. 
7 


P(X,r) =F(X) +r 
i=l j=l 


The general solution procedure is then to seek minimum values of P(X,r) 
money — 1, ,t5,...,t, Where r, >ro... ty - 0. Phen 


lim X, ee 
k-+ © 


and 


lim P(X, ,t, ) =F 
k-+ @ 


where F is the minimum value of the objective function of the mathema- 
tical programming problem at x, 

The conditions for the determination of a global minimum for the P 
function are generally the same as for the convex programming problem. 
The objective function, F(X), must be a convex function. Recall that a 


linear function is convex although not strictly convex. Also the 


ey 


constraint region must be convex (which is not true for this problem) 
and F ,G 72... G, » eee a must be continuous. 


Close inspection of the P function reveals that the term 


is a boundary repulsion term which keeps the minimum feasible solution 
of the P-function in the interior of the region defined by the inequality 
constraints. Naturally, the use of such a boundary repulsion term 
requires that solutions exist in the interior of the region defined by the 
inequality constraints. 

It is this feature of the Fiacco-McCormick technique which provides 
a significant improvement over other convex programming techniques. 
Other approaches, including Schmit's "method of alternate steps" require 
elaborate techniques to determine what action to take when the boundary 
is encountered during the minimization process. 

Another consequence of the boundary repulsion term which should be 
noted, however, is that the search for a feasible P-function minimum must 
be confined to the region defined by the inequality contraints. The 
reason for this restriction is apparent from the following example. Suppose 
that G(X) = x, 2 0. As x,-0 from the negative side IG, )]7~ -@. 

The P-function can be minimized by direct analytical procedures for 
relatively trivial cases only. For other cases, some method of successive 
approximations must be used. One such method is that of gradient 


descent. 
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A. GRADIENT DESCENT 
Gradient methods are based ona Taylor's series approximation of 
the function to be minimiZed. 
] — Eirst-Order Grddienmiethed 
The first-order gradient method uses only the first-order terms of 
the Taylor's series expansion to obtain the equation, 
Xp =X, - 8 7 P(X) 
The procedure is then to move a distance proportional to @ in the direc- 
temo tie gradient’vector evaluated at the point X,. x, 19 the pemmiat 
which the P-function is minimized along the gradient. 
2. Second-Order Gradient Method 
A better approximation of the P-function is obtained by using the 
second-order terms of the Taylor's series expansion as well as the first 


to obtain the equation 


=i 
a. | 7 P(X) 


== 4 = @ | 
OX, OX, 





where the matrix ||- || is the Hessian matrix evaluated at X, and -# 
denotes the matrix inverse. The P-function is minimized along the vector 
given by the product of the inverse of the Hessian and the gradient of the 
function evaluated at X,. 

If the function to be minimized is not strictly convex, the Hessian 
matrix may be positive definite only in the vicinity of a local minimum. 


For this reason, the procedure used began with a search in the direction 


of positive 8. If the P-function did not decrease in this direction before 


encountering the boundary, the search was continued in the direction of 


negative @. This situation actually occurred in the example when the 
P-function was minimiZed at r=l. 
oe) blgher-@rdcmiMetiods 

It would seem logical to consider the use of higher order terms to 
obtain an even better approximation of the function to be minimized. The 
effort required to program such a procedure, however, would be prohibitive 
for all but trivial problems. In fact, the labor required to calculate the 
second-order partial derivatives for the Hessian matrix may prohibit the 
use of the second-order gradient method for large problems. Experience 
with several small examples using both techniques indicates that the 
convergence properties of the second-order technique are far Superior to 
those of the first-order technique, Fiacco and McCormick provide exam- 
ples of both techniques in Ref. 3 which clearly demonstrate the super- 


iority of the second-order method. 


B, ESTIMATION OF THE FINAL P-FUNCTION MINIMUM 

Fiacco and McCormick show in Ref. 2 that estimates of the solution 
of the P-function minimum for r=0 can be obtained by the use of poly- 
nomial approximations of x,,Xz,...,x, and P as a function of rt)? wae 
follows. Expand each component of the X vector associated with the 
minimum of P(X,r) for a given value of r in a power series in r*’*» To 
obtain a (p-1) th-order polynomial approximation, drop the terms of the 


power series after the pth term. If the P-function has been minimized for 


r,, k=l,...,h where h 2p, form the linear equations: 


70 


x, _ =a. eee Ay (, ) +...+ 


x (r)) =a +yanit+ Aya (0,) aE Seen aets! “ ) 


If it is assumed that ris reduced by a constant factor such that 
eee = Vi YA SCN rt, =T, To", then the above equations constitute a system 


oe linear CGudtions in Pp UNnKkKmawns , "de,... 72 For the polynomial, 


mele 


p-1 


os 2 
xy =a Ae ar ayt aia ast ar ® rc) e AF eee: i 


aeeat r—0 is given by a. Hence, if these equations are Solved for am, 
a (p-1) th-order polynomial approximation is obtained for the value of x, 
corresponding to the solution to the mathematical programming problem. 
Eoumates are obtained by this procedure for x,,... 1X as well as P at 
r=0. Because of computer accuracy limitations, the order of the poly- 
Hoemial estimate will be restricted to three and will use the last tour 


data points. 


C. CRITERIA FOR TERMINATION 

As r approaches zero, the value of the P-function minimum approaches 
the optimal solution. More rapid convergence to the optimal solutiom may 
be obtained by the extrapolation described above thereby permitting term- 
ination of the procedure with fewer iterations. The specific point at 
which the minimiZation process is to terminate depends, of course, upon 


the accuracy desired of the final solution estimate. 


Zl 


Fiacco and McCormick suggest several alternate criteria in Ref. 2 
based on the changes in the values of the variables or functions from one 
P-minimum to the next. In the case of the estimated value of the P- 
function minimum, the minimiZation process is terminated when 

|P*(ry_.) - P¥ (1,) |< 
where € is a small positive number and P*(r,) is the polynomial estimate 
of the final solution after the P-function is minimized for the ith r value. 


For the three-bar truss problem, € was chosen to be 1xl10°°, 


D,. SUMMARY OF SOLUTION PROCEDURE 
The following procedures are those required to obtain an estimated 
value for the optimal solution to the mathematical programming problem. 

1. Form the P-function and obtain the analytical forms for the gradient 
vector and Hessian matrix, 

2. set p=4 for a third-order polynomial estimate of the values of P and 
Kye eK for r=0. 

3, oetr=l1 and h=1 and select ajsiarting point, ~j.such that thesia. 
equality constraints are satisfied. 

4. Determine the direction of steepest descent at the starting point. 
search along this direction in the region defined by the inequality con- 
Straints until the directional derivative of thesP=function is lesismtham the 
specified tolerance, i.e. near zero. 

9. Test the gradient magnitude. If less than a specified tolerance, go 


to step 8. 


Za 


6. Call the present point the new starting point and go to step 4. 

7. Ifh 2p calculate the (p-1) th-order polynomial estimate of the com- 
ponents of the solution at r=0. 

8, If h 2 p+l compare the last two estimates of the P function at r=0. 
If within the specified tolerance, terminate the routine. Otherwise, set 


r =r/C where C is the reduction factor and seth = h+l. Goto step 6. 


Bi) 


V. THE EXAMPLE 
The example used to illustrate the application of the Fiacco- 
McCormick technique to the design of trussed structures is the problem 
described previously with two distinct load conditions. This yields a 
problem with 13 variables; 10 equality constraints, 4 of which are non- 


linear; and 23 inequality constraints including the three non-negativity 


constraints for the cross-sectional areas of the truSs members. 


A, PROGRAMMING CONSIDERATIONS 

The computer program was written in Fortran IV language for the 
IBM 360 computer at the Naval Postgraduate School in Monterey, 
California. No attempt was made to optimize the program to reduce 
execution time. As it is written, the program required approximately 25 
seconds for execution. 

A simple bracketing procedure was used to locate the minimum of the 
P-function along the vector in the direction of steepest descent. The 
directional derivative was evaluated at the current starting point. A @é 
was then determined corresponding to a point at which the directional 
derivative was opposite in sign from that at the starting point and still 
within the region defined by the inequality constraints. The bracketing 
procedure was then commenced to reduce the magnitude of the directional 
derivative below the specified tolerance, DDTOL. 

Some experimentation was required to determine the necessary mag- 
nitude of DDTOL to ensure convergence of the gradient magnitude to within 


a tolerance of EPSI = 1xl10~°. This was determined to be DDTOL = 1xi0°’. 
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For the three-bar truss problem, a double-precision version of the 
program was required although several two-variable problems were 
solved without recourse to double precision. For small r values, the 
P-function becomes progressively more peaked in the vicinity of the 
minimum value. For the three~bar truss problem this peakedness is such 
that the change in @required to obtain a bracket small enough to yield a 
directional derivative within tolerance is below the accuracy of the 
computer unless double precision is used, 

For programming purposes, the notation of the problem variables 


WV A, ,Bs,As »0739 70a 951, Se. 97 2,055,904 0,U,,U 
gO BAB OLLI RL ISB LI Py te) 26 B23 219!" 2 ~was changed to 


Bec, se, &. .,X, 5 TeESpectively:. 


B. INPUT DATA 

The input data for the problem was that of Schmit's Case 1 in Ref. l. 
A cursory inspection will reveal that the values are highly unrealistic 
since the actual range of values of the modulus of elasticity is from 5 to 
30x10° lbs./in.* for metals. They are, however, sufficient for purposes 


of illustrating the solution technique. The input data is as follows. 


P, = 30 0, = 1 f, = 135° E, = 1 
pe oom be = ee E, =a 
@,= 60° 9, =1 B; = 45° =, =a 
a. = 180° 
Ves ol 


The constraints specified for the problem are 


ae, 


eee ae 
-15 = 01, ,%2;293 | = yA BS ile, Zug 
and 


=|5 (arse nuy 


tae = 200,e9=1..2" 
xj) Yj) 


The starting point for the example was arbitrarily chosen to be xi = 


Cele BUSeSeseoeoes enon s)— The reduction factor for r was 10 at 


each iteration, i.e. r4,, =r, /10. 


CORED ULES 

Values obtained for minimum P as a function of r as well as the third- 
order estimates of the P-function minimum for r=0 are shown in Table I. 
The value of the objective function, F, corresponding to the indicated r 
value; the total value of the remaining terms of the P function, denoted 
by P-F; and the values of the equality constraints, H,(X) are also shown, 

The value given under "Moves Required to Minimize" is the number 
of times a new direction of steepest descent was calculated prior to 
achieving a gradient magnitude of less than 001 for the current r value. 

The worth of the polynomial estimate is indicated by the fact that a 
final solution estimate of minimum P is obtained on the sixth iteration 
which differs by less than 1x10°* from the value obtained on the ninth 
and final iteration. 

By the ninth iteration, the value of the objective function, F, has 
converged to within 7x10-° of the final solution estimate while the value 


of P(r.) differs by 14x107° from the final estimate. 
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One consequence of not obtaining an absolutely null gradient vector 
is revealed by the values of the equality constraints. Since convergence 
is obtained only to a gradient magnitude of .001, the equality constraints 
are not exactly satisfied for large r values. The values of the equality 
constraints are seen to decrease, however, at each iteration. In other 
words, the constraints get progressively "tighter" as r~ 0. The reason 
for this is the peakedness of the P-function in the vicinity of the minimum 
as r gets smaller and smaller. As P becomes more peaked, it is necessary 
to get nearer the minimum to satisfy the gradient magnitude tolerance 
for convergence. Thus the approximations become more accurate as 
r > 0. A gradient magnitude tolerance of 1x107~* failed to yield a 
change in the value of the constraints of sufficient magnitude to change 
the results although it did require a slightly higher execution time. 

Babble show satm@eavclltcs Omtne SuCceoclVemat), 1albaauwerotcs 
well as the final solution estimates. The final solution estimates changed 
by less than 2x107~* from the sixth to the ninth iteration. 

The values obtained by the method of Fiacco and McCormick are 
compared with those obtained by Schmit [1] using his "method of alter- 


nate steps" in Table III. 


oil 


TABLE Til 


Fiacco & Alternate 
Variable McCormick steps 
A, Pore i. Owe 
A, ,044 .944 
iL 
A, 61a Pour 
O11 Poe O77 19.863 
oo) Zo. JOG 19.984 
O14 mAs: 1.206 
u 19.755 19.743 
ai 
u 20.000 -19.984 
yl 
O15 -15.000 -14,993 
O45 9.000 S. UUS 
O 
32 20.000 19.996 
u -35 .000 -34.998 
Boe 
u - 9.000 - o, 003 
y2 
W DEES Lie, ae 24 


The results obtained by the two methods are quite close although the 
minimum weight determined by the method of Fiacco and McCormick was 
Slightly less. The primary reason for the difference was probably the 
termination criterion. Schmit solved the problem in 1960 on an IBM 653 
digital computer with severe storage capacity limitations so the accuracy 
of the solution in this study is unquestionably higher. Advances in 
hardware make any comparison of execution times meaningless. Schmit's 
results for o3, is erroneous and is probably a typographical error in the 


reference, 
a 


VI. DISCUSSION 


A, ALTERNATE MINIMA 

A solution to the three-bar truss problem has been determined but 
the fact remains that it can be verified to be only a local minimum, Pro- 
cedures exist by which the "goodness" of the solution may be determined 
but the effort may be more than that required for solving the entire problem 
by a trial and error solution technique. One such procedure is simply to 
vary the starting point and search for alternate minima, 

The example problem was solved for several other starting points 
selected at random. These included 
Jt 


X = (110), 1.0,(y o- OO. won 


OR OIS 20M 


a (0.25,0.25,0.25,0.0,...,0.0)" 


= (0.25 ORZ SDE? Sean COR OO 
= (0025 ,0.25 07 Seta ee 


The only change in value noted was a difference in the number of 


moves required to minimize the P-function for r=l. 


B, TWO-BAR TRUSS 

If the center bar of the truss is removed, the structure becomes a 
statically determinant two-bar truss in which the internal forces are 
uniquely determined. The minimum weight of such a structure subject to 


the constraints of the example is 3.04905. Thus, by adding a third 
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member, the total weight of the structure is reduced. It should be noted, 
however, that the behavior of the trusses will not be the same, 

Switsky [4] shows that the deflection of the two-bar truss will be 
less in this case than that of the three-bar truss. If the displacements 
of both trusses are required to be equal, given the same stress con- 
straints, the statically determinant two-bar truss will be the lighter of 


the two. 


C. NON-ANALYTICAL CONSTRAINTS 

Structural problems often do not allow the nice problem statement 
form of the three=bar truss problem. The major difficulty arises from not 
being able to express some of the constraints in the completely analytical 
form required for the Fiacco~-McCormick technique, Many problems 
require the use of empirically determined alignment charts or graphs to 
calculate changes in structural behavior as a result of the variation of 
design variables. For such problems, some variation of Schmit's 
"method of alternate steps" may be well suited. An example of such a 
problem is the determination of the minimum weight design of an integrally 
stiffened panel loaded axially under compression. Rosenbaum [5] des- 
Ccribes the problem and outlines a solution procedure based on the method 


of alternate steps, 
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VII, RECOMMENDATIONS FOR FURTHER STUDY 

Since only three design variables of the three-bar truss are allowed 
to vary, the example is actually a sub-optimization problem. This fact 
Suggests that several possible extensions of the problem merit further 
consideration, 

Allowing more design variables such as the length of truss members 
to vary makes the problem even more nonconvex but the Fiacco-McCormick 
procedure may Still be applied, The assumption of the availability of a 
continuous spectrum of materials would permit the determination of 
optimal material properties of the structure. 

The added restriction of a discrete choice of materials or material 
sizes available is a further realistic extension of the problem although a 
different solution procedure is required, 

Finally, the development of a systematic procedure for finding alter- 


nate minima is essential if one is to accept any result with confidence. 
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